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We construct time reversal invariant topological superconductors and superfluids in two and 
three dimensions which are analogous to the recently discovered quantum spin Hall and three-d Z2 
topological insulators respectively These states have a full pairing gap in the bulk, gap less counter- 
propagating Majorana states at the boundary, and a pair of Majorana zero modes associated with 
each vortex. We show that the time reversal symmetry naturally emerges as a supersymmetry, 
which changes the parity of the fermion number associated with each time-reversal invariant vortex. 
In the presence of external T-breaking fields, non-local topological correlation is established among 
these fields, which is an experimentally observable manifestation of the emergent supersymmetry. 
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The search for topological states of quantum mat- 
ter has become an active and exciting pursuit in con- 
densed matter physics. The quantum Hall (QH) effect [l| 
provides the first example of a topologically non-trivial 
state of matter, where the quantized Hall conductance 
is a topological invariant Q . Recently, the quantum spin 
Hall (QSH) stated, i 

has been theoretically predicted 
Q and experimentally observed in HgTe quantum well 
systems [6|. The time reversal invariant (TRI) QSH state 
is characterized by a bulk gap, a Z 2 topological number 
0, and gapless helical edge states, where time-reversed 
partners counter-propagate [jj 0]. 

Chiral superconductors in a time reversal symmetry 
breaking (TRB) (p x -\-ip y ) pairing state in 2d have a sharp 
topological distinction between the strong and weak pair- 
ing regimes [Io|. In the weak pairing regime, the system 
has a full bulk gap and gapless chiral Majorana states at 
the edge, which are topologically protected. Moreover, a 
Majorana zero mode is trapped in each vortex core[10j, 
which leads to a ground state degeneracy of 2 n_1 in the 
presence of 2n vortices. When the vortices wind around 
each other a non-Abelian Berry phase is generated in the 
2 n_1 dimensional ground state manifold, which implies 
non-Abelian statistics for the vortices [11]. Chiral super- 
conductors are analogous to the QH state — they both 
break time reversal (TR) and have chiral edge states with 
linear dispersion. However, the edge states of a chiral su- 
perconductor have only half the degrees of freedom com- 
pared to the QH state, since the negative energy quasi- 
particle operators on the edge of a chiral superconductor 
describe the same excitations as the positive energy ones. 

Given the analogy between the chiral superconducting 
state and the QH state, and with the recent discovery of 
the TRI QSH state, it is natural to generalize the chiral 
pairing state to the helical pairing state, where fermions 
with up spins are paired in the (p x + ip y ) state, while 
fermions with down spins are paired in the (p x — ip y ) 
state. Such a TRI state have a full gap in the bulk, 
and counter-propagating helical Majorana states at the 
edge (in contrast, the edge states of the TRI topologi- 
cal insulator are helical Dirac fermions). Just as in the 



case of the QSH state, a mass term for an odd number 
of pairs of helical Majorana states is forbidden by TR 
symmetry, and therefore, a topologically protected super- 
conducting or superfluid state can exist in the presence 
of time-reversal symmetry. Recently, a Z2 classification 
of the topological superconductor has been discussed in 
Refs [12, LL3|, LL4J, by noting the similarity between the 
Bogoliubov-de Gennes (BdG) superconductor Hamilto- 
nian and the QSH insulator Hamiltonian. The four types 
of topological states of matter discussed here are sum- 
marized in Fig. [TJ In this work, we give a Z2 classifica- 
tion of both the 2D and 3D cases which has a profound 
physical implication. In two dimensions, we show that 
a time-reversal invariant topological defect of a Z2 non- 
trivial superconductor carries a Kramers' pair of Majo- 
rana fermions. Let Np be the operator which measures 
the number of fermions of a general system, then the 
fermion- number parity operator is given by (— 1)^ . This 
operator is also referred to as the Witten index [151]. which 
plays a crucial role in super symmetric theories. We prove 
the remarkable fact that in the presence of a topological 
defect, the TR operator T changes the fermion number 
parity, T~ 1 (-1) Nf T = -(-l) Np locally around the de- 
fect in the Z2 non- trivial state, while it preserves the 
fermion number parity, T _1 (— 1) Nf T = (—1) Nf , in the 
Z2 trivial state. This fact gives a precise definition of the 
Z2 topological classification of any TRI superconductor 
state and is generally valid in the presence of interactions 
and disorder. A supersymmetric operation can be de- 
fined as an operation which changes the fermion number 
parity; therefore, in this precise sense, we show that the 
TR symmetry emerges as a supersymmetry in topolog- 
ical superconductors. Though supersymmetry has been 
studied extensively in high energy physics, it has not yet 
been observed in Nature. Our proposal offers the op- 
portunity to experimentally observe supersymmetry in 
condensed matter systems without any fine tuning of mi- 
croscopic parameters. The physical consequences of such 
a supersymmetry is also studied. 

As the starting point, we consider a TRI p-wave su- 
perconductor with spin triplet pairing, which has the fol- 
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FIG. 1: (Top row) Schematic comparison of 2d chiral super- 
conductor and the QH state. In both systems, TR symmetry 
is broken and the edge states carry a definite chirality. (Bot- 
tom row) Schematic comparison of 2d TRI topological super- 
conductor and the QSH insulator. Both systems preserve TR 
symmetry and have a helical pair of edge states, where oppo- 
site spin states counter-propagate. The dashed lines show 
that the edge states of the superconductors are Majorana 
fermions so that the E < part of the quasi-particle spectra 
are redundant. In terms of the edge state degrees of freedom, 
we have (QSH) = (QH) 2 = (Helical SC) 2 = (Chiral SC) 4 . 



lowing 4x4 BdG Hamiltonian: 
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with #(x) = ^C|(x), C|(x), c|(x), cj(#)^ , e p = p 2 /2m — 
\± the kinetic energy and chemical potential terms and 
h.c. = ( i (72 cr a A a ipjY. The TR transformation is de- 
fined as Cf —> C|, ci — > — C|. It can be shown that the 
Hamiltonian (pQ) is time-reversal invariant if A a j is a real 
matrix. To show the existence of a topological state, con- 
sider the TRI mean-field ansatz A al = A(l, 0, 0), A a2 = 
A(0, 1,0). For such an ansatz the Hamiltonian ([1]) is 
block diagonal with only equal spin pairing: 
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with S[f(x) = ^C|(x), cj(x), C|(x), c|(x)^ , and p± = 

± From this Hamiltonian we see that the spin 
up (down) electrons form p x + ip y (p x — ip y ) Cooper 
pairs, respectively. In the weak pairing phase with 
H > 0, the (p x + ip y ) chiral superconductor is known 
to have chiral Majorana edge states propagating on 
each boundary, described by the Hamiltonian if e dge = 
^2k y >o v Fky^-k y ^k y ^eie ip- ky is the quasipar- 

ticle creation operator [10[ and the boundary is taken to 
be parallel to the y direction. Thus we know that the 
edge states of the TRI system described by Hamiltonian 
([2]) consist of spin up and spin down quasi-particles with 



opposite chirality: 

Hedge = Vpk V (^"M^M - ^-kyi^kyi) ■ (3) 



k v >0 



The quasi-particle operators ^k y h V^l can be expressed 
in terms of the eigenstates of the BdG Hamiltonian as 

^k y \ = J d 2 x (u ky {x)c^{x) + v ky (x)c\(x)^ 

^k y [ = J d 2 x (u*_ ky (x) Ci (x) +v*_ ky (x)c\(x)^j (4) 

from which the time-reversal transformation of the quasi- 
particle operators can be determined to be T~ 1 ?p ky ^T '■ 



il>-k v i, T 1 ^k y [T -i/j-kyP In other words, 

(ipkyh ^-kyi) transforms as a Kramers' doublet, which 
forbids a gap in the edge states due to mixing of the 
spin-up and spin-down modes when TR is preserved. To 
see this explicitly, notice that the only /^-independent 
term that can be added to the edge Hamiltonian (J3]) is 
im^2 k ip-ky^ipkyi with m G R. However, such a term is 
odd under TR, which implies that any back scattering be- 
tween the quasi-particles is forbidden by TR symmetry. 
The discussion above is exactly parallel to the Z2 topo- 
logical characterization of the quantum spin Hall system. 
In fact, the Hamiltonian ([2]) has exactly the same form as 
the four band effective Hamiltonian proposed in Ref. Q to 
describe HgTe quantum wells with the QSH effect. The 
edge states of the QSH insulators consist of an odd num- 
ber of Kramers' pairs, which remain gapless under any 
small TR-invariant perturbation [gl, A no-go theorem 
states that such a "helical liquid" with an odd number 
of Kramers' pairs at the Fermi energy can not be real- 
ized in any bulk Id system, but can only appear as an 
edge theory of a 2d QSH insulator [9[. Similarly, the edge 
state theory ([3|) can be called a "helical Majorana liquid" , 
which can only exist on the boundary of a Z2 topolog- 
ical superconductor. Once such a topological phase is 
established, it is robust under any TRI perturbations. 

The Hamiltonian (pQ) can be easily generalized to three 
dimensions, in which case A aj becomes a 3 x 3 matrix 
with a = 1,2,3 and j = x,y,z. An example of such a 
Hamiltonian is given by the well-known 3 He BW phase, 
for which the order parameter A aj is determined by an 
orthogonal matrix A aj = Au a ^ u G SO(3)[l6j]. Here and 
below we ignore the dipole-dipole interaction term [ItJ 
since it does not affect any essential topological proper- 
ties. By applying a spin rotation, A aj can be diagonal- 
ized to A aj = A5 a i , in which case the Hamiltonian (pQ) 
has the same form as a 3d Dirac Hamiltonian with mo- 
mentum dependent mass e(p) = p 2 /2m — \i. We know 
that a band insulator described by the Dirac Hamilto- 



nian is a 3d Z2 topological insulator for \i > 18.ll9l.l2 



and has nontrivial surface states. The corresponding su- 
perconductor Hamiltonian describes a topological super- 
conductor with 2d gapless Majorana surface states. The 
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surface theory can be written as 

#surf = \ ^2 ^F^-k (<Tzk x +&xk y ) ^k (5) 
Z k 

which remains gapless under any small TRI perturba- 
tion since the only available mass term m ^-k^y 
is time-reversal odd. We would like to mention that the 
surface Andreev bound states in 3 He-B phase have been 
observed experimentally [2lj . 

To understand the physical consequences of the non- 
trivial topology we study the TRI topological defects of 
the topological superconductors. We start by consid- 
ering the equal-spin pairing system with BdG Hamil- 
tonian ([2]) in which spin up and down electrons form 
Px + i>Py and p x — ip y Cooper pairs, respectively. A 
TRI topological defect can be defined as a vortex of 
spin-up superfluid coexisting with an anti- vortex of spin- 
down superfluid at the same position. In the generic 
Hamiltonian (pQ), such a vortex configuration is written 
as A aj = [exp(ia 2 <9(r-r ))] aj ', a = 1,2 and A 3 ^' = 0, 
where 0(r — tq) is the angle of r with respect to the vor- 
tex position ro. Since in the vortex core of a weak pair- 
ing p x + ip y superconductor there is a single Major ana 
zero mode[10|,l22|, one immediately knows that a pair a 
Major ana zero modes exist in the vortex core we study 
here. In terms of the electron operators, the two Majo- 
rana fermion operators can be written as 

7, - /A(«(x)c,(x) + „;(x)c| W ) 

7t =/A(«5(x)c t W + » (x)c !W ) (6) 

where we have used the fact that the spin-down zero 
mode wave function can be obtained from the time- 
reversal transformation of the spin- up one. The Ma- 
jorana operators satisfy the anti-commutation relation 
{lailp} — 28 a p. The TR transformation of the Majo- 
rana fermions is 



-7 T . 



(7) 



Similar to the case of the edge states studied earlier, the 
Major ana zero modes are robust under any small TRI 
perturbation, since the only possible term imj^Ji which 
can lift the zero modes to finite energy is TR odd, i.e., 

The properties of such a topological defect appear iden- 
tical to that of a 7r-flux tube threading into a TRI topo- 
logical insulator 23l 24|, where a Kramers' pair of com- 
plex fermions are trapped by the flux tube. However, 
there is an essential difference. From the two Majorana 
zero modes 7t , 7j a complex fermion operator can be de- 
fined as a = (7t +*7|)/2, which satisfies the fermion 
anticommutation relation {a, a''"} = 1. Since 7^,7^ are 
zero modes, we obtain [a,H] = which implies that a 



is the annihilation operator of a zero-energy quasipar- 
ticle. Consequently, the ground state of the system is 
at least two- fold degenerate, with two states \Gq) and 
\Gi) — o)\ Gq) containing and 1 a- fermions. Since 
a^a — (1 + Z7|7|) /2, the states |G?o(i)) are eigenstates of 
^7t7l w ^h eigenvalues — 1(+1), respectively. Thus from 
the oddness of ij^Ji under TR we know that \Go) and 
\d) are time-reversed partners. Note that superconduc- 
tivity breaks the charge £7(1) symmetry to Z2, meaning 
that the fermion number parity operator (—1) Nf is con- 
served. Thus, all the eigenstates of the Hamiltonian can 
be classified by the value of (—1) Nf . If, say, \Gq) is a 
state with (—1) Nf = 1, then \Gi) = \Gq) must sat- 
isfy (—1) Nf = —1. Since \Go) and \Gi) are time-reversal 
partners, we know that in the Hilbert space of the zero- 
energy states the TR transformation changes the fermion 
number parity: 



t-\-i) Nf t = -(-i 



(8) 



Eq. dSJ is the central result of this paper. At a first 
glance it seems contradict the fundamental fact that the 
electron number of the whole system is invariant under 
TR. Such a paradox is resolved by noticing that there are 
always an even number of topological defects in a closed 
system without boundary. Under the TR transforma- 
tion, the fermion number parity around each vortex core 
is odd, but the total fermion number parity remains even 
as expected. Once the anomalous transformation prop- 
erty (jSJ) is established for a topological defect in a TRI 
superconductor, it is robust under any TRI perturbation 
as long as the bulk quasiparticle gap remains finite and 
other topological defects are far away. Thus Eq. (|8j) is a 
generic definition of TRI topological superconductors: 

• Definition I. A two-dimensional TRI supercon- 
ductor is Z2 nontrivial if and only if fermion num- 
ber parity around a TRI topological defect is odd 
under TR. 

A transformation changing fermion number by an odd 
number is a "super symmetry" ; thus, the TR symmetry 
emerges as a discrete supersymmetry for each TRI topo- 
logical defect. The same analysis applies to the edge 
theory (|3]), which shows that in the Id helical Majorana 
liquid is a theory with TR symmetry as a discrete super- 
symmetry. 

All the conclusions above can be generalized to 3d 
topological superconductors. In the 3 He BW phase the 
Goldstone manifold of the order parameter is A aj = 
Au aj e £0(3) x 17(1) [S SH. A time-reversal invari- 
ant configuration satisfies A aj G R, which restricts the 
order parameter to SO (3). Since Hi(SO(3)) = Z2, the 
TRI topological defects are Id "vortex" rings. By solv- 
ing the BdG equations in the presence of such vortex 
rings, it can be shown that there are linearly dispersing 
quasiparticles propagating on each vortex ring, similar 
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to the edge states of the 2d topological superconductor. 
However, for a ring with finite length the quasi-particle 
spectrum is discrete. Specifically, there may or may not 
be a pair of Majorana modes at exactly zero energy. The 
existence of the Majorana zero modes on the vortex rings 
turns out to be a topological property determined by the 
linking number between different vortex rings. Due to 
the length constraints of the present paper, we will write 
our conclusion and leave the details for a separate work: 
There are a pair of Majorana fermion zero modes con- 
fined on a vortex ring if and only if the ring is linked to 
an odd number of other vortex rings. Such a condition 
is shown in Fig. [2l Consequently, the generalization of 
Definition I to 3d is: 

• Definition II. A 3d TRI superconductor is Z 2 
nontrivial if and only if the fermion number parity 
around one of the two mutually-linked TRI vortex 
rings is odd under TR. 



(a) 



(a) 



(b) 




FIG. 2: Illustration of a 3d TRI topological superconductor 
with two TRI vortex rings which are (a) linked or (b) un- 
linked. The E — k\\ dispersion relations show schematically 
the quasiparticle levels confined on the red vortex ring in both 
cases, "o" and "x" stand for the quasiparticle levels that are 
Kramers' partners of each other. Only case (a) has a pair of 
Majorana zero modes located on each vortex ring. 

Besides providing a generic definition of the Z 2 topo- 
logical superconductors, such an emergent super symme- 
try also leads to physical predictions. Consider the 
2d topological superconductor coupled to a weak TR- 
breaking field M(r), which is classical but can have ther- 
mal fluctuations. This situation can be realized in an 
isolated superconductor with vortices pinned to quenched 
weak magnetic impurities. The n-point correlation func- 
tion of M(r) can be obtained by 
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(9) 



in which the trace is restricted to states with an even 
number of fermions. For a closed system with N vor- 
tices, the leading order effect of the TR-breaking field is 
to lift the degeneracy between the two Majorana fermions 
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FIG. 3: (a) Illustration of a 2d TRI topological supercon- 
ductor with four TRI topological defects coupled to a TR- 
breaking field. The arrows show the sign of the TR-breaking 
field sgn(M(r s )) at each topological defect. In the two con- 
figurations shown, only the field around vortex 1 is flipped, 
leading to an opposite fermion number parity in the corre- 
sponding ground state \G) and \G') (see text), (b) Illustra- 
tion showing the flow of the energy levels when the upper 
configuration in figure (a) is deformed to the lower one. The 
flip of the TR-breaking field M(ri) leads to a level crossing at 
M(ri) — 0, where the fermion number parity in the ground 
state changes sign. 



in each vortex core. Consequently, the perturbed Hamil- 
tonian H[M(r)] to first order can be written as 



H[M(r)} 



N 



iM(r s )a s 7474 



(10) 



in which J s m) are the- Majorana fermion operators, and 
a s G K depend on the details of the perturbation. The 
important fact is that the mass term induced is linear in 
M(r) at the defect position r s , since ijs^jsi is TR odd. 

Since the superconductor has a full gap, naively one 
would expect all the correlations of M(r) field to be short 
ranged. However, for a system with N topological defects 
the TV-point correlation function has a long range order 
when r s , s = 1, 2, ..N are chosen to be the coordinates of 
the topological defects. In other words, the correlation 
function 



, f im v . Vf[M(r s )Wo, 

|ri-rj|-K», V*,j \ " / 



(ii) 



though all the n point correlations in Eq. © with n < N 
remain short ranged. Physically, such a non-local correla- 
tion can be understood by comparing two states \G) and 
|G'), which are the ground states of the systems with the 
field configurations M = (M(ri), M(r 2 ), M(r N )) and 
M' = (-M(ri),M(r 2 ), ...,M(r N )), respectively. From 
Hamiltonian ([TO)) it can be seen that \G) and \G r ) have 
opposite fermion number parity, since the fermion num- 
ber parity around the first topological defect Z7i|7ij is 
reversed while that of all the other topological defects 
remains invariant, as shown in Fig. [3l Without loss of 
generality, we can assume ( — 1)^ is even for \G) and odd 
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for \G'). Since the whole system is required to have an 
even number of fermions, the lowest energy state in the 
Hilbert space for the field configuration M 1 is not |G'), 
but the lowest quasiparticle excitation c^^^G'). Thus, 
the two field configurations M. and M 1 have different 
free energies, which leads to the non- vanishing correla- 
tion function in Eq. ([TT]) . Even when the topological 
defects are arbitrarily far away, the energy difference be- 
tween the two configurations remains finite, which shows 
the non-local topological correlation. Similar non-local 
correlations can also be obtained for a 3d TRI topologi- 
cal superconductor with linked vortex rings. 
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